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Abstract 

Through a finite size renormalization group technique we calculate 
the running coupling constant for quenched SU(2) with a few percent 
error over a range of energy varying by a factor thirty. The definition is 
based on ratio of correlations of Polyakov loops with twisted boundary 
conditions. The extrapolation to the continuum limit is governed by 
corrections due to lattice artifacts which are proportional to the square 
of the lattice spacing and appears rather smooth. 
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1 Introduction 

Lattice gauge theories can challenge the determination of the strong coupling 
constant obtained from high energy experiments. The most precise data are 
from the LEP collaborations who derive a value for a s with a six percent 
error at the Z energy scale. This error turns into a fourty percent error on 
A-qcd because asymptotic freedom makes low energy values of the coupling 
constant very sensitive to their high energy values. Precise tests of the stan- 
dard model require a higher precision for a s , of the order of one percent at 
the Z mass. Indeed, global fits to LEP data show a variation of the strong 
coupling constant with the top mass of the order of one to two percent for a 
twenty GeV variation of the heavy quark mass. In this paper we present a 
lattice gauge theory calculation of the running coupling constant in quenched 
SU(2) which shows that the required precision can indeed be reached. The 
method follows the work of refs. [1, 2, 3, 4] where a finite size scaling tech- 
nique was first applied to the estimate of the running coupling constant. 
Our work differs from the ones quoted above for the choice of the observ- 
able which defines the renormalized coupling. The definition of refs. [2, 3, 4] 
is based on the Schrodinger functional while ours is related to the ratio of 
correlations of Polyakov loops with twisted boundary conditions. The choice 
of twisted boundary conditions was dictated by the necessity of performing 
a weak coupling (small g 2 ) lattice perturbative calculation to determine the 
relation between our definition and the other ones. Indeed, according to the 
definition, the calculation must be performed on a finite size lattice where 
twisted boundary conditions eliminate the problem of torons. This calcula- 
tion was completed and the results reported in ref. [7]. Section 2 contains the 
generalities of the definition and reviews the method. Section 3 discusses the 
details of the simulation and the use of lattices with an odd number of points 
in spite of the fact that our definition refers to correlations of Polyakov loops 
taken at a distance of half the lattice size. Section 4 contains the results for 
the running coupling constant including the extrapolation to the continuum 
limit and the calibration of the energy scale. A first successful comparison 
of these results with those of the Schrodinger functional is made. A detailed 
comparison will be made elsewhere [8]. Finally, we present the results for the 
/3-function. 



2 The definition and the method 

The calculation adopts the following boundary conditions: 

U^r + uL) = 0^(r)fit v = x, y 

U^r + uL) = U^r) v = z,t (1) 

with U^(r + uL) a "replica-link", i.e. a link in a site displaced by the lattice 
size in the u direction and Q x , Q y two unitary matrices satisfying appropriate 
commutation rules (see refs. [5, 6] for details). We have imposed twisted 
boundary conditions on two directions which eliminate the finite volume 
degeneracy of toron configurations with the perturbative vacuum and allow 
a standard perturbative expansion [7]. Polyakov loops which wrap around 
the lattice in the twisted directions (e.g. the a;-direction) are defined by: 

L-l 

P x (y,z,t) = Tr([n U x (x=j,y,z,t)]Sl x e i 5v) (2 ) 

3=0 

where Q x is the twisting matrix in the direction of the loop. This form is de- 
termined by the requirement of gauge and translation invariance. The above 
definition produces correlations of "twisted" loops which start at order a s in 
the weak coupling limit. In addition, the disconnected term in the correlation 
is identically zero. However, such a correlation cannot be taken straightly 
as a definition of the coupling because of linear divergences, characteristic 
of Wilson lines, which spoil the continuum limit. A proper definition is ob- 
tained by dividing such a correlation by a similar one between Polyakov loops 
extended in a direction with periodic boundary conditions (z or t). 

A crucial ingredient of the method is the definition of the energy scale 
of the renormalized coupling related to the total lattice size, which allows 
a finite size scaling recursive scheme. The distance between Polyakov loops 
should be taken equal to some fixed fraction of the lattice size, half in our 
case. By defining the coefficient of the tree level result extrapolated to the 
finite-volume continuum limit as T(oo), the definition for g 2 (aL) reads: 

2, n 1 (E y , z P x (0,0,0)PKy,z,t = L/2)) 

q (aL) = — — 1 (3) 

V ^ T(oo) (£,,, P z (0, 0, 0)P2(x, y , t = L/2)) V ' 

where the correlation in the denominator has to be understood as a full one, 
disconnected part included. 



The ratio of the A parameter of this definition (Atp) to the MS one, 
calculated in ref. [7], is given by: 

A rP /A S 5 = 1.6136(2) (4) 

By replacing the normalization with the tree level result T(oo) at L = oo 
in eq.(3) with the tree level T(L) at finite Z, one obtains a definition which 
differs by terms of order a 2 and therefore is characterized by the same Ktp- 
The recursive scheme, already explained at length in the original papers, 
will be only summarized here: 

• on a volume with L points per side, and physical linear size L p h ys = La 
at a value of (3 = 4/ ' g\ large enough to be close to the continuum limit, 
one calculates the renormalized coupling g 2 (La); 

• by repeating the same calculation for the case with 2L points, one 
gets the coupling g 2 {2La) renormalized at a length scale which is twice 
larger; 

• the previous two steps are repeated for larger values of (3, i.e. for 
smaller values of the lattice spacing, while the total initial physical 
volume Lp h is kept constant. This is achieved by changing the number 
of lattice points and the value of the bare coupling so that the value of 
the renormalized coupling stays constant; 

• the set of calculations at different lattice spacings are extrapolated to 
zero lattice spacing giving the value of the change of the renormalized 
coupling by a change in the length scale by a factor 2 in the continuum; 

• the continuum value of the renormalized coupling at the larger scale is 
taken as the input value for a new step and the procedure is iterated 
until one reaches a physical volume which is large enough to make the 
calculation of a physical quantity, which will calibrate the value of the 
lattice spacing, free from finite size effects; 

• the normalization of the lattice spacing on the largest volume sets the 
scale of the running coupling constant in physical units. 



3 The details 

The simulation was made by using an action with ordinary periodic boundary 
conditions and a tower of frustrated plaquettes in the x-y plane, which can 
be shown [6] to be equivalent to a standard action with the twisted boundary 
conditions defined in eq.(l). 

The updating scheme (HFKP) was the modification of the standard Creutz's 
heat bath algorithm introduced in refs. [9, 10] which improuves the accep- 
tance of the original algorithm at high values of (3. In order to fight the 
critical slowing down at large volumes and high (3, the HFKP algorithm was 
followed by a number of overrelaxed iterations increasing with the volume 
up to 12 for the case L = 16 at (3 = 3.08 . The Polyakov loops measurement 
is affected by large fluctuations: the number of measurements, the auto- 
correlation time t and the relative error e obtained for different values of 
(3 and different volumes are given in table 1. In order to reduce the 

fluctuations, we have used the "one link integral" technique by replacing the 
product of links in eq. (2) with the product of the values of the links averaged 
in the "heat bath" of their neighbours. For the SU(2) group one obtains the 
following analytic expression: 

itt t \\ kl {k) - 2h{k) — t 

mx))o LI = ^ V,(x) (5) 

where k = f3 detfT^z)) 1 / 2 , V^x) = V^(x)/ det(V^(x))^\ V^x) is the "sta- 
ple" relative to the link U^(x) and I n is the modified Bessel function of order 
n. The variance of the correlation of the twisted Polyakov loops improuves 
with the "one link integral" technique by a factor which increases from 2 to 
3 for (3 decreasing from 4 to 3. 

In spite of the definition which involves the time distance L/2, we did also 
use lattices with an odd number of points. Consider our definition of g 2 (aL) 
(eq. (3)) calculated at time distance (L — l)/2 in lattice units. Because 
of periodic boundary conditions in the time direction such a definition is 
actually the arithmetic average of the one at distance t = (L — l)/2 and 
of the one at distance t = (L + l)/2 and therefore differs from the one at 
t = L/2 by an even function of 1/ L. This guarantees that it will have the 
same A and then the same continuum limit. However, in order to plot the 
results of odd lattices on the same extrapolation curve of even lattices, we 
need a definition which reproduces the one at t = L/2 up to terms at least 
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Table 1: The number of over-relaxation sweeps per heat-bath sweeps N 
the number of measurements N meas , the auto-correlation times r and tl 
relative errors e obtained for different values of (3 and L 



of order 1/L 4 . This is possible if our definition admits a Taylor expansion 
around t = L/2. Both numerical and perturbative investigations show that 
this is the case. A linear combination can then be formed of operators at 
t = (L — l)/2 and t = (L — 3)/2 which has the desired property that its 
expectation value differs from the one of the operator at t = L/2 by terms of 
order 1/L 4 only. If we rewrite the definition at t = L/2 (eq. 3) in the form 

2, n 1 h<fi,L,L/2) 

9 ( ai > = rR /«(/*,£, 1/2) ' < 6 » 

where the expectation values f x and / z , which depend on the bare coupling 
constant /3, the lattice size L and the correlation time distance L/2, are de- 
fined by comparison with eq. (3), the required linear combination is obtained 
by making in eq. (6) the following replacements: 

/.(/?,£,f) - 9 -f x ((3,L,^)-^f x ((3,L,^) 

/,(/?, £,f) - 9 -f z ((3,L,^)-^f z ((3,L,^) (7) 

We will see that points obtained with this definition do align with those for 
even lattices in the curves showing the extrapolation to zero lattice spacing 
limit and indirectly confirm the smallness of order 1/L 4 corrections. 



4 The extrapolation to the continuum limit 

Getting rid of the lattice cutoff is the most delicate aspect of the calculation. 
Table 2 contains the results at fixed renormalized coupling for the value of 
the coupling renormalized at a double length scale, for different values of the 
lattice spacing, i.e. on a fixed physical volume for different values of L. The 
same results are reported in figure 1. The extrapolation to the continuum 
limit appears rather smooth. We outline the points obtained on odd lattices 
which align with the even lattice points. The value of the coupling at the 
scale Lphys should be kept constant in the extrapolation procedure, but this 
is possible only within the statistical errors. In order to obtain the value 
of (3 which corresponds to the reference value of g 2 (L p h ys ) we have made an 
interpolation among the (3 values nearby and the resulting error is taken as 
the uncertainty on the value of the reference coupling. By using a two loop 



p 


L 


g\La) 


g 2 (2La) 


3.9762 


4 


2.075 


2.334(12) 


4.0300 


5 


2.075 


2.381(11) 


4.0873 


6 


2.075 


2.410(15) 


4.1780 


8 


2.075 


2.417(22) 


3.7109 


4 


2.430 


2.848(10) 


3.7620 


5 


2.430 


2.942(18) 


3.8210 


6 


2.430 


2.942(14) 


3.9200 


8 


2.430 


2.975(16) 


3.3946 


4 


3.068 


3.809(14) 


3.4563 


5 


3.068 


3.853(23) 


3.5179 


6 


3.068 


3.892(19) 


3.6200 


8 


3.068 


3.975(34) 


3.0879 


4 


4.108 


5.527(22) 


3.1607 


5 


4.108 


5.552(37) 


3.2235 


6 


4.108 


5.719(27) 


3.2770 


7 


4.108 


5.602(49) 


3.3300 


8 


4.108 


5.621(48) 


2.8349 


4 


5.597 


8.065(42) 


2.9103 


5 


5.597 


8.126(61) 


2.9745 


6 


5.597 


8.161(48) 


3.0305 


7 


5.597 


8.207(69) 


3.0800 


8 


5.597 


8.039(87) 



Table 2: The renormalized couplings for different (3 and L 
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Figure 1: The extrapolation to the continuum limit 
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Table 3: The values for g 2 {2L p h ys ) compared with the two-loops ones 

expression for the /3-function, the errors on the coupling at scale L p h ys are 
converted into errors on the coupling at the scale 2L p h ys and are added to its 
statistical fluctuations as additional and independent statistical fluctuations: 
in table 2 the errors on the coupling at the scale 2L p h ys include these two 
contributions. 

Table 3 shows the numerical values extrapolated to the continuum limit 
of g 2 (2L p h ys ) and a comparison with the values predicted by the two-loop 
evolution for that step. The continuum limit is obtained by making a least- 
squares fit of the data for g 2 {2La) with a curve linear in L~ 2 . 

Figure 2 shows the values of the renormalized coupling extrapolated to 
the continuum limit as a function of an energy scale in arbitrary units. 

Besides the error on the extrapolation to the continuum limit of table 3, 
there is a second source of errors when the different values of the renormalized 
couplings of table 3 are put in sequence to reconstruct the running coupling 
constant. Indeed, each recursive step has an error on the coupling at the large 
scale which must be propagated through the various steps. This is done by 
converting the uncertainty on the coupling g 2 {2L p h ys ) into an uncertainty 
on the renormalization scale which then would fluctuate around an exact 
multiple of two. However, the renormalized coupling at the scale L p h ys of the 
step n-\- 1 is not exactly equal to the one at the scale 2L p h ys of the step n. By 
using the two-loops expression for the /3-function, this small difference can 
be turned into a slight modification of the renormalization scale which will 
fluctuate at the various steps around a value different from an exact multiple 
of two. 

This is reported in table 4, by taking as the reference value for the scale 
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Figure 2: The running coupling constant as a function of the energy scale in 
arbitrary units 
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Table 4: The running coupling constant as a function of the scale given in 
units of the reference scale L and in GeV 

the one of the fifth iteration step, where we will perform the calibration of 
the scale in physical units. In figures 2 and 4, the errors on the scale are 
turned again into errors on the couplings. These figures contain points at 
scale lengths which are without error while a cumulative error bar appears 
on the values of the coupling. The dashed line in figure 2 is the behaviour 
predicted by the two-loops continuum /3-function normalized to the coupling 
at highest energy which closely follows the running coupling up to the last 
points and justifies the use of the two-loops /3-function in the error propa- 
gation procedure described above. 

The calibration of the length scale, i.e. of the lattice spacing, requires 
a precise calculation of a physical quantity on the largest volume and close 
enough to the continuum limit. In our case we have not performed such 
a dedicated calculation which is certainly crucial for the real unquenched 
SU(3) case. Following ref. [3] we have taken the available values of the string 
tension at a (3 value 2.85 overlapping with the values of (3 of the fifth step 
of the recursive procedure. This can be seen from figure 3 which shows the 
dependence of (3 upon L at fixed renormalized coupling: the line shows a 
fit with the dependence implied by asymptotic scaling which is in very good 
agreement with the results of the simulation. The calibration should be taken 
as an exercise rather than as a precise determination: indeed such a (3 value 
lies between the calculations of the renormalized coupling on a L = 4 and 
an L = 5 lattice. The mild L dependence of our results for the renormalized 
coupling justifies the utility of such an exercise. A careful normalization will 
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Figure 3: The dependence of (3 upon the number of points L at constant 
renormalized coupling g 2 (L p h ys ) = 5.597 
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Figure 4: The running coupling constants denned in this paper (o) and the 
one of ref. [3] (□) as a function of the energy scale with the same lattice 
calibration 



be discussed elsewhere [8]. By adopting for the string tension a conventional 
value of (425 MeV) 2 we get that the scale for the running coupling constant 
shown in figure 4 and in the third column of table 4 is in GeV. The lower 
points in figure 4 correspond to the calculation performed in ref. [3] with the 
same scale calibration. The two definitions should match in the perturbative 
region, once after implementing the relative correction term. This can be 
included by rescaling the calibration of one scale with respect to the other 
by the ratio of the A parameters given by: 



Atp/Asf = 3.6206(2) 



(8) 
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Figure 5: The running coupling constants denned in this paper (o) and the 
one of ref. [3] (□) with the same lattice calibration after rescaling the energy 
scale by the ratio of A parameters 

where Asf is the A parameter of the definition based on the Schrodinger 
functional (see ref. [2]). 

Figure 5 shows the result of this rescaling: at small values of the coupling 
the two definitions nicely superimpose to each other and, partly surprisingly, 
they do so also at rather large values of the coupling. 

A result independent from A is the /3-function: by locally fitting the 
running coupling constant with a polynomial one can extract the /3-function 
for our definition and the one of ref. [3]. As an exercise we fitted three points 
at time with a second-order polynomial in \og(E): the result is plotted in 
figure 6, together with the two-loops /3-function. Beyond a value of the 
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Figure 6: /3-function for the running coupling constants denned in this paper 
(o) and the one of ref. [3] (□) is compared with the two-loops expression 
(dotted line) 

coupling around 6 our /3-function deviates from the two-loops result. 

5 Conclusions 

The above simulation shows the feasibility of a calculation of the running 
strong coupling constant at a few percent level. The calculation took about 
1,600 CPU hours at an average speed of about two Gigaflops on a three Gi- 
gaflops version of the APE100 series. The unquenched case at least requires, 
with the present algorithms, a factor hundred in computer time. Moreover, it 
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requires the understanding of the inclusion of fermions with twisted boundary 
conditions for which interesting suggestions have been made. 
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